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ABSTRACT 

We present deep HST/NICMOS Camera 3 F110W and F160W imaging of a 26'x33', corresponding to 3.1pcx3.8pc, non-contiguous 
field towards the Orion Nebula Cluster (ONC). The main aim is to determine the ratio of low-mass stars to brown dwarfs for the cluster 
as a function of radius out to a radial distance of 1.5pc. The sensitivity of the data outside the nebulous central region is F160W< 21.0 
mag, significantly deeper than previous studies of the region over a comparable area. We create an extinction limited sample and 
determine the ratio of low-mass stars (0.08-1 M Q ) to brown dwarfs (0.02-0.08M o and 0.03-0.08M o ) for the cluster as a whole and 
for several annuli. The ratio found for the cluster within a radius of 1.5pc is R 01 = N(0.08 - lM e )/N(0.02 - 0.08M o ) = 1.7 ± 0.2, and 
R 03 = N(0.0S - 1M )/A'(0.03 - 0.08M o ) = 2.4 ± 0.2, after correcting for field stars. The ratio for the central 0.3pcx0.3pc region down 
to 0.03M o was previously found to be R 03 = 3. 3^ 7, suggesting the low-mass content of the cluster is mass segregated. We discuss 
the implications of a gradient in the ratio of stars to brown dwarfs in the ONC in the context of previous measurements of the cluster 
and for other nearby star forming regions. We further discuss the current evidence for variations in the low-mass IMF and primordial 
mass segregation. 

Key words, open clusters and associations: individual (ONC) stars: formation stars: low-mass, brown dwarfs stars: luminosity 
function, mass function stars: pre-main-sequence 



1. Introduction 

A central question in astrophysics is whether the Initial Mass 
Function (IMF) is universal or if it is a function of star formation 
environment. The answer is crucial for our understanding of the 
buildup of galaxies and provides clues to how stars and brown 
dwarfs are formed. The stellar part of the IMF a ppears to be uni- 
versal based on the comparison of the field I MF dChabrierll2003l) 
to clus ters within 2 kpc (e.g. iMever et ail |2000|) . iBastian et all 
d2010l) reached the same conc lusion with the possi ble excep- 
tion of the Taurus Dark Could (iLuhman et al.L I2009). For more 
massive Galactic clusters and for low-metallicity clusters in 
the Large Magellanic Clouds the IMF is well represented by 
a po wer-law to the sensitivity and crowding limit of ~1M^ 
(e.g.lStolte et all l2Q06h iGouliermis et al.Ll2007tlBrandner et all 
2008: 1 Andersen et alll2009l) . 

Th e situation is l ess cl ear fo r the brown dw arf part of the 
IMF. iBriceno et al. | d2002l) and iLuhmanl d2004l) suggested a 
higher ratio of stars to brown dw arfs in Taurus than the Orion 
Nebula Cluster (ONC). However. iGuieu et all ri2006l) utilized a 
larger area survey of Taurus and found the ratio o f stars to brown 
dwarfs to be similar for th e two regions (c.f. ILuhman et all 
2009i). lAndersen et alJ {2008) combined the measured IMF's for 
six nearby star-forming regions and the Pleiades and showed 
that there is no strong evidence for variations in the derived ratio 
of low-mass stars to brown dwarfs for the regions and that the 
combined ratio for the 7 regions showed that the underlying IMF 
was falling into the brown dwarf regime. Conversely, a handful 
of studies for OB associations over large areas of the sky have 
indicated an IMF that does not show eviden ce for a falling mas s 
function in the brown dwarf regime (e.g. lLodieu et all 20071) . 



The reason for the apparent discrepancy is not clear but could 
be due to either differences in the star formation process or due 
to dynamical evolu tion of the clusters. For several of the clus- 
ters compiled by Ander sen et alJ ((2008) only the central regions 
were studied (e.g. a 0.6pc square box for the ONC and 0.24pc 
square box for Mon R2) and the measured IMF in the centre 
can thus be biased if mass segregation is present. Large field of 
view surveys of young star forming regions are needed to di- 
rectly compare with the OB associations. 

One nearby region of particular interest is the ONC. So 
far, the cluster content has only systematically been stud- 
ied dee p into the brown dwarf regim e for the central ~5' 
square (Iffillen brand & Carpenterl l2000t iMuench et ail I2002L 
ISlesnick et all |2004|) . corresponding to ~ 0.6pc square for 
a distance of 414pc (|Menten et al., 2007). Larger area, less 
deep, studies of the stellar population have indicated at least 
the higher mass objects in th e ONC is mass segregated 
(Hill enbrand & HartmannLll998l) . Thus, the low-mass IMF de- 
rived for the central area can be biased towards higher mass ob- 
jects and may not be representative of the cluster as a whole. 
However, no studies have combined deep observations in the 
near-infrared together with a large field of view coverage such 
that the brown dwarf content can be studied in detail as a 
function of radius and to determine the global ratio of stars 
to brown dwarfs. Here we present HST NICMOS Camera 3 
F110W and F160W (similar to the ground-based J and H bands) 
non-contiguous observations of the ONC area. The observations 
cover the cluster from the very centre out to a radius of 2.55pc 
and were obtained as part of the HST Orion Treasury Project 
(GO 10246, PI Robberto). 
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The paper is structured as follows. Section 2 presents the 
observations and data reduction. The source detection and the 
photometry is discussed in Section 3. In Section 4 we derive the 
ratio of low mass stars to brown dwarfs as a function of radius in 
the cluster and we discuss the evidence for the lower mass clus- 
ter content to be mass segregated. In Section 5 we compare the 
findings for the ONC with the similar results from other nearby 
star forming regions. Finally, we conclude in Section 6. 



2. Observations 

Regions of the ONC were observed using the HST/NICMOS 
Camera 3 during the HST Treasury Program covering the Orion 
Nebula Cluster (GO10246). The non-contiguous nature of the 
NICMOS Camera 3 observations stems from the fact they were 
observed in parallel with the ACS and WFPC2 Treasury Project 
observations. The field-of-view of the NICMOS Camera 3 is 
51'.'2x51'.'2 and the pixel scale is 072. Since historically, the 
NICMOS Camera 3 is slightly out of focus, there is a slight 
degradation of the point spread function (PSF) which in turn im- 
proves the sampling of the PSF. 

A total of 102 regions were covered in both the F110W and 
F160W bands. The regions were typically build up as a 4 piece 
mosaic in the F160W band and a 5 piece mosaic in the F110W 
band utilizing the Field Offset Mechanism (FOM, see below). 
In one point, a single 5 172x5 172 image was observed in both 
filters. The integration time of each image was 192 seconds in 
F160W and 256 seconds in F110W, resulting in a total integra- 
tion time per pixel that varies between 192-762 seconds for the 
F160W observations and 256-1280 seconds for the F110W ob- 
servations, with the central parts of each sub-mosaic being im- 
aged for the full amount of time. Around 55% of the mosaics are 
covered for 192 seconds in the F160W band, and 40% for 380 
seconds. For the F110W filter, around 45% is covered for 256 
seconds, and 40% for 512 seconds. 

The large dithering pattern was possible due to the FOM, 
which however introduces some vignetting. Fortuitously, a 3x3 
NIC MOS Camera 3 mos aic of the Trapezium was obtained in 
1998 (lLuhmanetal.Ll2000t) . We were able to make ratios of one 
vignetted FOM-dithered exposure to that unvignetted exposure. 
A smooth function was fit to the vignetted edge in the ratio image 
and then expanded into a 2-D vignetting correction image (one 
per filter for F110W and F160W) that was applied to all affected 
exposures. The raw images have further been processed with the 
appropriate dark and flat frames as described on the STScIs web 
pages. They were further corrected for bad pixels determined 
from sigma clipping of the median stacked image. 

Offsets were determined between the individual images of 
each region by identifying stars in common and determining the 
centroid. The frames were then combined using the drizzle soft- 
ware that takes the distortion of the NICMOS Camera 3 frames 
into account and combines the frames with sub-pixel accuracy. 
All the source detection, astrometry, and photometry is being 
done on the final mosaic for each region. 



3. Source detection and photometry 

We describe the source detection and photometr y process. The 
photom etry is compared to the previous work by Rob berto et alJ 
(1201 Oh and is converted into the 2MASS system. Completeness 
corrections are determined and the fractional coverage of the 
ONC is described. 



3.1. Source detection and photometry 

Source detection is complicated by strong nebulosity, especially 
close to the centre of the ONC, and by the under-sampled PSF of 
the NICMOS Camera 3. A ring median filtered image was sub- 
tracted from each mosaiced image in order to remove the large 
scale nebulosity. The inner and outer annuli for the ring filter 
were 7 and 11 pixels, (174 - 272). This makes source detection 
easier but the photometry was performed on the frames with the 
nebulosity maintained. The PSF is better sampled in the F160W 
images than the F110W images due to the larger full width at 
half maximum at longer wavelengths. Since the F160W data are 
also deeper than the F110W data, we have used the F160W data 
for the source detection. 

Point sources were visually identified in each individual 
frame. The radial profile of each source was inspected in order 
to exclude spurious detections due to nebulosity. We detect a to- 
tal of 2054 objects. The positions of the objects detected in the 
F160W images were then used for the photometry in the F110W 
band. The photometry was performed with an aperture radius 
of 2.5 pixels (075) and a sky annulus between 10 and 15 pixels 
(2"-3") for both filters. Zero points were adopted from the HST 
calibration web pages. Bright isolated stars in different regions 
were used to determine aperture corrections. The random error 
on the aperture correction was 0.02 mag for both bands. The 
adopted zero points, including the aperture corrections, are 22.50 
mag for the F110W band observations and 21.66 mag for the 
F1 60W band obse r vation s. Objects in common with the dataset 
by iRobberto et al.l d2010h were used to calibrate the astrometry 
of each mosaic. For one of the mosaics, no ISPI star was present 
within the field-of-view. Coordinates for the stars in this field 
were obtained from the header information and is expected to 
be accurate to ~ 1", the average accuracy of the header coor- 
dinates for the mosaics that could be calibrated with ISPI. The 
astrometry is in agreement with the ISPI astrometry to within 
071 for 60% of the sources and within 0/2 for almost 90% of 
the sources. The sources and their photometry are presented in 
TableQ] 

3.2. Integrity of the photometry 

The NICMOS photometry has been calibrated into the 2MASS 
system utilizing the ISPI data which was transformed into the 
2MASS system dRobberto et all 1201 Oh . We chose relatively well 
exposed sources, not located close to the edge of the NICMOS 
mosaics. Least square linear fits were performed to the F160W 
band magnitudes and the Fl 10W-F160W color. 

H 2 mass = F160W + a 0H (1) 
Jimass = F110W + aoj + a u * (FU0W - F160W), (2) 

where aon = -0.25 + 0.01, any = -0.26 + 0.01, and au = 
-0.35 + 0.01, respectively. Fig. Q] shows the ISPI photometry 
as a function of the NICMOS photometry for both the J and H 
bands. All sources with a nominal photometric error smaller than 
0.1 mag in each band both in this study and the ISPI observa- 
tions are shown. Further, the photometric error for the NICMOS 
sources is shown for both filters. The agreement with the ISPI 
photometry is in general good but there is some scatter. This is 
to be expected due to the difference in spatial resolution making 
it easier to separate starlight and nebulosity in our observations 
but also due to our poorly sampled PSF, especially in the / band. 
In any case, our photometry is significantly deeper than the pre- 
vious ISPI study, which reached 3cr limits of J~ 19.5 mag, and 
H~ 18 mag. The faintest sources identified in this study have 
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Table 1. The sources detected in the NICMOS Camera 3 mosaics. Objects located in frames that could not be astrometrically 
calibrated are marked. Magnitudes are in VEGAMAG. The full table is electronically available on A&As web site. 



ID 


RA 


DEC 


F110W 


eFHOW 


F160W 


eF160W 




(J2000) 


(J2000) 


(mag) 


(mag) 


(mag) 


(mag) 





83.55908 


-5.49581 


23.71 


0.80 


21.23 


0.36 


1 


83.55791 


-5.49974 


22.53 


0.44 


19.84 


0.18 


2 


83.55621 


-5.50233 


19.18 


0.09 


17.09 


0.05 


3 


83.55345 


-5.50198 


23.53 


0.70 


21.12 


0.33 


4 


83.55041 


-5.50002 


23.13 


0.58 


20.45 


0.24 


5 


83.54137 


-5.50148 


22.86 


0.53 


20.62 


0.27 


6 


83.54234 


-5.49987 


23.26 


0.65 


20.50 


0.25 


7 


83.54798 


-5.50549 


14.19 


0.01 


13.06 


0.01 


8 


83.54616 


-5.50673 


22.85 


0.52 


21.60 


0.43 


9 


83.54810 


-5.50938 


20.40 


0.16 


18.29 


0.09 



0.35-0-7 pc 



0.7-1-1 pc 






20 25 





Fig. 1. Upper and lower left: A comparison of the NICMOS pho- 
tometry converted into the 2MASS system wit h the photome- 
try ob tained with ISPI, in the 2MASS system (Rob berto et al.L 
2010). Only objects with photometric errors less than 0.1 mag 
in both filters are shown. Upper and lower right: The NICMOS 
photometric errors as a function of magnitude. 



Fig. 2. The completeness correction for different annuli as a 
function of F160W band magnitude. From the upper left is 
shown the completeness corrections for the inner annulus (0.35- 
0.7pc), the middle annulus (0.7-l.lpc), the outer annulus (1.1— 
1.5pc), and in the lower-right a comparison of the three annuli. 



magnitudes of J~ 22.0 mag, and H~ 21.0 mag for the areas 
with limited nebulosity. The discussion below describes in more 
detail how the completeness limit depends on the distance from 
the cluster centre, and hence on the amount of nebulosity. 

3.3. Completeness corrections 

We have performed completeness experiments to determine the 
fraction of stars detected as a function of object magnitude. 
Bright isolated stars from several frames were used to create 
a model for the point spread function. Artificial point sources 
were then placed randomly within the observed frames, with 
F160W magnitudes randomly chosen between 15 and 20 mag. 
In each frame around 80% of the detected number of stars were 
added as artificial stars. The frames are not crowded and the ar- 
tificial star experiments are being performed to test the detection 
limit as a function of nebulosity. The relatively large number of 
stars added per frame does not influence this limit. A color of 
Fl 10W-F160W-2 mag, was adopted for all the artificial stars, 
corresponding to the colour of the faintest objects within the ex- 
tinction limited sample defined below (see Fig . O ■ 

The artificial stars were recovered by hand in the F160W 
frames in a similar manner to the real star detection. Photometry 
was then performed in both the F160W and the F110W frames. 



An artificial star was considered detected if the magnitude error 
in both filters was less than 0.3 mag, similar to the maximum 
error for some objects used in the extinction limited sample. 

Fig. |2] shows the completeness corrections for three differ- 
ent annuli as a function of F160W band magnitude. The 50% 
completeness limit clearly varies across the survey area. The 
presence of strong nebulosity close to the cluster centre results 
in a brighter 50% completeness limit, as expected. In particu- 
lar, within the central 0.35pc radius, the strong nebulosity and 
undersampled PSF makes the detection of fainter point sources 
difficult and the completeness corrections are substantial. On the 
other hand, farther than 0.35pc from the centre the complete- 
ness is high down to F160W-12, mag. For comparison, a 1 Myr 
20Mj„p object at the distance of the ONC has H band magnitude 
of 16.4 mag. Thus, even at 0.35pc we can detect essentially all 1 
Myr objects down to 20My„ p seen through a n extinction of Ay = 
7 mag (see below), adopting the models by Bara ffe et alj <[l998). 
At larger radial distances objects below 10M /U/J are detected. 
Due to the incompleteness within 0.35pc we focus our analysis 
on the objects at larger ra dii. The current survey thus com ple- 
ments the obs e rvatio ns by Hillen brand & Carpenter! (2000) and 
iLuhman et al.1 (120001) that both cover the central 0.6pc square 
region. For all radii the completeness limit is higher than the de- 
tection limit of the data, which is due to both the nebulosity and 
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0.35-0.70 pc 



35m36.0s 12.0s 48.0s 24.0s 5h34m00.0s 



Fig. 3. F160W band 26'x33' mosaic of the observed region 
within the ON C. For refer e nce, th e central 140"xl40" region 
observed by iLuhman et al.l d2000l) is included. The location of 
the Trapezium stars are marked by small boxes in the mosaic. 
North is up and East is to the left. 




Distance from center (pc) 



Fig. 4. The fraction coverage of the ONC by the NICMOS 
Camera 3 observations. The fraction covered is shown as a func- 
tion of radial distance from the centre {8 l Ori C). 
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Fig. 5. The color-magnitude diagr am for the surveyed region. 
Overplotted are a 1 Myrisochrone (Baraff e et al.L 1 1 998b shifted 
to the distance of the ONC, the same isochrone reddened by 
Ay =7 mag, and the 50% completeness limit for the data. 
Masses for a 1 Myr population of 1.0, 0.08, 0.02, and 0.01M o 
are indicated on the lower left panel. The region used as control 
field is shown in the lower right corner (2.3-2.6pc). The typi- 
cal photometric and colour error are shown in each diagram as a 
function of object magnitude. 



4. Results and Analysis 

We now present the J-H versus H color-magnitude diagrams for 
the surveyed region and estimate the contamination due to field 
stars. The mass of each object is derived from its de-reddened 
H band magnitude. We create an extinction limited sample to 
ensure an unbiased selection of cluster members as a function 
of object mass before determining the ratio of low-mass stars to 
brown dwarfs for the ONC as a function of radius. 



to the structure in the nebulosity that can hide in particular faint 
sources. 

3.4. Fraction of the ONC covered 

Since the observations were obtained in parallel, the field-of- 
view covered is not contiguous. The very centre of the cluster 
is not covered and the majority of the observations are to the 
south and south-south-west of the cluster centre. Fig. [3] shows 
the location in the ONC of e ach individual mosaic . For refer- 
ence, the mosaic observed by Luhmar Tet alJ (|2000) is inserted 
and the Trapezium stars marked. Fig.|4]shows the fractional area 
of the ONC covered as a function of distance from the cluster 
centre. 

The total area imaged is 162 square arcminutes, correspond- 
ing to 2.35 pc 2 for a cluster distance of 414pc. The coverage 
extends in radius from 0.05pc out to 2.55pc from 9 l Ori C. Due 
to the large region of sky covered with a small filling factor, the 
coverage is relatively low, 20-50% within 0.5pc, 10-30% be- 
tween 0.5pc and 1.5pc, and less at larger distances as shown in 
Figg] 



4.1. The color-magnitude diagram 

Fig. [5] shows the color-magnitude diagram for the surveyed re- 
gion outside 0.35pc. The data have been divided into the annuli 
that will be used below to derive the ratio of low-mass stars to 
brown dwarfs. Overplot ted are the 50% co mpleteness limits, a 
1 Myr isochrone from iBaraffe et al.l d 19981) co nverted into the 
2MAS S system using the same procedu re as in | Andersen et al.1 
(2006) with the filter transformations by ICarpent er (2001), and 
the same isochrone reddened by Ay = 7 mag. 

Most of the objects from 1M to the detection limit below 
0.0 1M are located between the 1 Myr isochrone shifted by the 
distance modulus to the ONC, and the same isochrone reddened 
by Ay = 7 mag. The objects on the blue side of the 1 Myr 
isochrone are likely foreground stars or a significantly older pop- 
ulation. 

4.2. Extinction limited samples 

One complication in calculating the IMF in young star forming 
regions is the presence of differential extinction. Objects are red- 
dened by different amounts due to patchy obscuration along the 
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Fig. 6. The derived extinction for sources in the ONC within the 
magnitude limits 10.5 < H < 18 mag, which corresponds to 
the range from where the ratio of stars to brown dwarfs is deter- 
mined. 



line of sight. As evident from Fig. [5] this is the case for the ONC. 
Higher mass objects can be detected through larger amounts of 
extinction than lower mass objects due to their i ntrinsic bright- 
ness. As illustrated in e.g. lAndersen et al.1 d2009t) this can po- 
tentially introduce a spurious deficit in the derived number of 
low mass objects. One approach to avoid this bias is to construct 
an extinction-limited sample, which is obtained by considering 
only objects within a certain range of extinction while taking the 
completeness limit into account. 

The lower mass limit for the extinction-limited sample is set 
by a combination of depth of the data, completeness corrections, 
and maximum extinction allowed. The choice of an extinction 
limit is a compromise between completeness and the low mass 
limit (lower extinction) and maximizing the sample size (higher 
extinction). 

The extinction for each individual object is determined by 
shifting the object in the color-magnitude diagr am back to the 
assumed isochrone, in a similar manner as in Andersen et al.l 
(2006). We find that more than 80% of the objects are affected 
by an extinction of Ay = 7 mag or less and use this limit for our 
extinction limited sample. 

We have examined whether the extinction distribution is sim- 
ilar at different distances from the centre of the ONC. Fig. [6] 
shows the cumulative distribution of derived extinction for ob- 
jects in the magnitude range 10.5 < H < 18 mag, similar to the 
magnitude range from where the IMF is derived. There is no ob- 
vious trend of the extinction distribution as a function of radius. 
A two-sided Kolmogorov-Smirnov test between the cumulative 
distributions for the different radii indicates that we cannot rule 
out that the underlying distribution is the same. The maximum 
deviations between the extinction distributions for different radii 
were 0.12 with the significance levels varying between 0.37 and 
0.62. Thus, since the distribution of extinction appears to be sim- 
ilar as a function of radius, the same extinction limit is used for 
all radii. 

The mass of each object is then determined by sliding it 
along the reddening vector in the colour-magnitude diagram 
back to the isochrone. Since only J and H band observations are 
available we have no knowledge if an object is associated with 
a circumstellar disc that would give rise to an infrared excess. 
With JHK photometry the disk is typically taken into account 
by de-reddenin g the objects to the T-Tauri locus in the colour- 
colour diagram (Me veret al.Lll997h 



Not taking disks explicitly into account will tend to over- 
estimate the derived mass compared to the actual mass. Thus, 
a strong radial dependence on the disk fraction could in prin- 
ciple mimic a decr easing ratio of st a rs to brown dwarfs as a 
function of radius. Hillenbran d et al. I (fl998h find a modest de- 
crease of the disk fraction, from some 65% showing strong evi- 
dence of a disk in the centre to 45% at a radius of l.lpc. The 
ratio determined from the cen t ral 0.3 pc suffers less from the 
effects of disks. ISlesnick et al.l ((2004) determined the spectral 
types through spectroscopy and although veiling is still a po- 
tential issue, the use of J band spectra reduces this effect since 
the photosphere for late type objects peak in this wavelength 
range. However, the bolometric magnitudes still depend on the 
de-reddening in the colour-magnitude diagram. Using different 
de-reddening approaches provide difference in the luminosity 
of AlogL = 0.05 and is not expected to have an influence on the 
ratio of stars to brown dwarfs. The colour-magnitude diagrams 
in Fig. [5] further indicate this is not likely to be the source of 
the decrease in the ratio of stars to brown dwarfs as a function of 
radius. We would expect to see the strongest effect is in the 0.35- 
0.7pc bin where few objects are located close to the star-brown 
dwarf division line. For objects at the star- brown dwarf bound - 
ary the excess is J-H of 0.1 mag or less ( Mey er et al J, 1 1997b . 
which after de-reddening will correspond to maximum shift in 
H band magnitude of 0.16 mag. At most 4 objects would be 
redefined from stars to brown dwarfs from 0.35-0.7pc assum- 
ing all objects had disks and all objects were affected by 0.16 
mag of excess. Several objects at larger distances could also po- 
tentially be re-classified, 1 object from 0.7-1. lpc, and 3 objects 
from l.l-1.5pc . Neither will change t he ratios significantly and 
according to lHillenbrand et all (1 19981) half or less of the objects 
would show evidence for an infrared excess. We thus conclude 
it is unlikely the decrease in the ratio of stars to brown dwarfs is 
due to a lower disk fraction at larger radii. 

4.3. Field star contamination 

Due to the relatively large area covered we expect some field 
star contamination. The contamination can be separated into two 
types, namely that of foreground and background stars. Since the 
ONC is located at only 414pc, we expect the foreground field 
star contamination to be relatively minor. In addition, the fore- 
ground objects will have relatively blue colors compared to the 
cluster members (see below). 

Background stars are a larger concern. Two effects limit 
the contamination. First, the ONC is located at (I, b) - 
(209.07, -19.49), out of the Galactic plane and away from the 
inner Galaxy. Further, the ONC and surrounding region is asso- 
ciated with and partly in front of a giant molecular cloud com- 
plex acting as a shield for b ackground objects . By comparison 
with the n CO line map from Ball v et alJ (1 19871) it is evident that 
all the observed region is located within the J shaped filament in 
Orion A-North, in agreement with strong emission seen in the 
IRAS 100 fim images. Based on the strength of the CO emission 
the extinction averaged over a 40"beam is well above Ay = 15 
mag, in agreement with the low number of reddened stars seen 
even far from the centre of the ONC. 

We have estimated the field contamination in the following 
way . Since the cluster surface dens ity at larger radii decreases as 
r" 2 dHillenbrand & Hartmanniri998h . we expect the outer region 
(2.0-2.6pc) to have a relatively small number of cluster mem- 
bers. Thus, the outer regions can be used as a conservative over- 
estimate of the field star contamination. However, as we indi- 
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cate below there appears to be an increase in the relative number 
of brown dwarfs at larger radii. Any cluster population present 
at large radii will therefore result in an over subtraction of the 
brown dwarf content. The derived ratio of stars to brown dwarfs 
after correction based on the outer regions will thus be an upper 
limit. 

A Galactic model could in principle be used to predict the 
field star contamin ation along the line of sight of the ONC (e.g. 
Robin et al., 2003). However, in practice this is difficult due to 
the lack of calibration of the models for faint magnitudes. The 
model can be used as a test of whether the control field re- 
produces the expected number of f oregro und field stars. Using 
the Galactic model by iRobin et all (120031) we find the predicted 
number of foreground field objects towards the control area (2.0- 
2.6pc) to be 8.9. The measured number of blue objects (J-H 
<0.8 mag) within the magnitude range 10 < H < 19 mag is 10, 
suggesting that the outer annulus traces the field star population 
well. Based on the agreement with the foreground population 
and the similarity in the extinction distribution as a function of 
radius, we adopt the outer annulus, 2.0-2.6pc as the control re- 
gion for both foreground and background stars. 

4.4. The low-mass Initial Mass Function 

From the extinction limited samples we can construct an unbi- 
ased estimate of the IMF down to the completeness limit of the 
data. As illustrated in Fig.[2]the completeness in our inner bin de- 
creases sharply below H=18mag but individual objects are still 
detected down to //=20mag. Thus, by restricting the sample to 
an extinction limit of Ay = 7 mag and a lower mass limit of 
20Mj up the sample is expected to be complete across the whole 
mass range for all radii. The completeness tests showed the in- 
ner annulus is incomplete for the 10-20Mj llp mass range and this 
part of the mass function is thus an upper limit. In addition, as 
described below, the field contamination is substantial for the 
lowest mass range. We thus mainly focus the discussion of the 
IMF down to 30Mj up but present the corresponding ratios for the 
samples exten ding to 10 an d 20Mj up . 

Following iMever et all (Hp00) and lAndersen etaD d2008l) 
we have determined the ratio of low-mass stars (O.O8-1M ) to 
brown dwarfs both for the cluster as a whole, and for 3 separate 
annuli outside a 0.35 pc radius. The results are summarized in 
Table [2] and the ratio of stars to brown dwarfs as a function of 
radius is shown in Fig. [7] where the ratios both with and without 
field stars subtracted are shown. 

The r atio determined at 0.3 5-0.7 pc is consistent with the 
results bv lSlesnicketail (2004) for the central 0.3 pc radius of 
the cluster. However, for larger distances there is a drop in the 
ratio of stars to brown dwarfs. Outside a radius 0.7 pc the ratio 
drops by almost a factor of two, although the error bars for the 
individual annuli are relatively large. The total ratio of stars to 
brown dwarfs for the cluster as a whole has been calculated by 
correcting the source counts in the different annuli by the frac- 
tion of the clus ter covered by that annulus (der ived from Fig. |4j). 
As shown by iHillenbrand & Hartmannl (Il998h the surface den- 
sity is well represented by a King profile, and thus extends as 
r~ 2 outside the core radius. Thus, even though the surface den- 
sity is lower at larger radii, this is compensated by an increase in 
the area in each annulus. The low ratio far from the centre drives 
the global ratio to a relatively low value of Rq^ = 2.4 ± 0.2, lower 
than determined in the cluster centre. 

The drop in the ratio of stars to brown dwarfs as the dis- 
tance from the centre increases is also seen for the sample with 
a mass limit of 20Mj up or l0Mj up instead of 30Mj up but due 
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Fig. 7. The ratio of low-mass stars (0.08-1. 0M Q ) to brown 
dwarfs (O.O3-O.O8M ) for the ONC as a function of radius. The 
plus signs with error bars indicate the field star corrected mea- 
surements whereas the squares indicate the measu rements before 
field s tar subtraction. The inner data point is from lSlesnick et ai] 
(2004) and the horizontal line is the average ratio with the lcr 
error interval provided by the dashed lines. 



to the stronger contamination of field stars as the sample is ex- 
tended to lower masses we focus the analysis on the higher mass 
limit data. Whereas the field contamination is some 10% for the 
brown dwarf content down to a mass limit of 30Mj up , it increases 
at lower masses and hence fainter magnitudes. This is a conse- 
quence of the underlying luminosity function for the field star 
population which extends as a power-law to faint magnitudes. 

Could the lower ratio of stars to brown dwarfs as a function 
of distance from the ONC centre be due to a ge effects? Ado pting 
the new lower distance to the ONC (Ment en et all 120071) and 
a more extensive sample of members whe re the spectral type 
has been determined through spectroscopy (Hillenbrand, 1997) 
or through narrow band imaging of TiO absorption features, 
iDa Rio etalJ (l2010h suggest that the cluster might be older than 
the previous estimate of 1 Myr. Recent work has further sug - 
gested an age gradient across the ONC (Reggiani et ail 120111) . 
where the centre is dominated by a younger population whereas 
a wider age distribution is seen at larger radii with older stars 
having a larger average distance from the centre. 

The effect of an older cluster at the same distance is that the 
pre-main sequence stars and brown dwarfs gets dimmer, an ef- 
fect that is stronger for lower masses and it could thus replicate a 
decrease of the ratio of stars to brown dw arfs if the age adopted 
for the outer parts of the cluster is too low. Reggi ani et al.l d201 ll) 
find an average age for the (stellar) population between 1 pc 
and 1.5 pc to be \og(age = 6.3). We have therefore also cal- 
culated the ratios of stars to brown dwarfs in the outer annuli for 
an age of 2 Myr. We find the ratio of low-mass stars to brown 
dwarfs only changes slightly from the values determined for a 1 
Myr isochrone, R m (2Myr,d = 1.1 - l.5pc)) = 1.1 + 08., and 
Ro 2 (2Myr,d = 1.1 - 1.5pc) = 0.8 + 0.6 and this will thus not 
change the conclusion that the ratio of stars to brown dwarfs de- 
crease as a function of radius. 



5. Discussion 

We present the constraints that our observations place on the 
global IMF for the ONC and compare the results with the 
field IMF and results from other nearby star forming regions. 
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Evidence for variations in the IMF as a function of environment 
is discussed. 

5.1. Variations in the ratio of stars to brown dwarfs as a 
function of radius 

The decreasing ratio of stars to brown dwarfs as a function of ra- 
dius suggests that the cluster is mass segregated. A sim ilar con- 
clusion was found by Hillen brand & Hartmannl d 19981) for the 
massive stars in the cluster centre. However, there was no in- 
dication for mass segregation for the lower mass population as 
their study only extended to O.3M and thus compared high mass 
stars to low-mass stars. In contrast, we find the sub-stellar mem- 
bers to be segregated compared to the stars as a whole, which 
are dominated by the lower mass members. 

iHillenbrand & Hartmannl d 19981) concluded that it is possible 
for the high-mass segregation to have an dynamical origin. This 
is not the case for the brown dwarfs observed here with the cur- 
rent cluster surfa ce density profile. The half mass relaxation time 
was estimated in Hillen brand & Hartmannl (11998) to be 6.5 Myr 
which is much longer than the age of the cluster (1-3 Myr) and 
thus unlikely to cause the variations in the ratio of stars to brown 
dwarf s through dynamical 2 body evolution (c.f. lBonnell et al.L 
119981) . 

A decrease of the ratio of stars to brown dwarfs as a function 
of radius has been observed in some open c lusters as well, e.g. 
IC 2391 dBoudreault & Bailer- Jonei 120091) where the brown 
dwarfs were found to have a flatter distribution than the low- 
mass stars. Similar l y, for IC 348 and central parts of the ONC, 
iKuma r & Schmeia (2007) use the minimum spanning tree to 
suggest the brown dwarfs are less centrally concentrated than 
the stars. A flatter distribution for brown dwarfs compared to 
stars is what would be expected from an ejection scenario where 
a fraction of the brown dwarfs are due to ejections in multi- 
ple systems during the earliest phases of star formation (e.g. 
iReipurth & daAeTl200Tl iBate et all 120031 

Mo re recent simulation by e.g. iBatd d2009t) : iBonnell et alj 
(1201 ll) does not predict an overabundance of brown dwarfs fur- 
ther away from the centre when the simulations end. However, 
the simulations typically end after 1-1.5 free fall times which, 
depending on the specific simulation corresponds to 0.2-0.5 Myr, 
significantly younger than the current ONC. Thus, a more ex- 
tended brown dwarf distribution could be generated by subse- 
quent dynamical interactions in the core leading to a substantial 
unbound population or a higher overall velocity dispersion. 

5.2. Comparison with other Star forming regions 

A smaller ratio of the number of stars to brown dwarfs for the 
whol e clus t er co mpared to the centre region has several implica- 
tions. lBatd (l2009l) predicted a value of R Q3 = 1 .25, based on SPH 
simulations of a spherical collapsing cloud producing a total of 
1254 objects, even lower than we determine here for the ONC. It 
is worth noticing that our value of /?03 = 2.4 + 0.2 is not consis- 
tent with a Chabrier IMF, which predicts too few brown dwarfs 
(R^iChabrier) = 4.95). The probability of the ratio to be drawn 
from an underlying population with a Chabrier distribution is 
less than 1%. The ratio we determine is in better agreement with 
a flat or rising IMF in linear units into the brown dwarf regime. 
We stress however, that our IMF is determined from a photomet- 
ric sample only and that spectroscopic confirmation of the brown 
dwarf candidates is highly desirable. 



lAndersen et alJ (120081) presented a compilation of well de- 
termined ratios of stars to brown dwarfs in nearby young star 
forming regions and the Pleiades. The measurements did not 
show evidence for a variable underlying IMF as a function of 
star formation environment. Results for some OB associations, 
e.g. Upper Sco, and cr Ori, have indicated that the derived mass 
function is r icher in brown dwarfs compared to the you ng cluster 
sample (e.g. lLodieu et al.L 120071; ICaballero et al.Ll2007l) . The de- 
rived low-mass slope is only weakly falling or flat into the brown 
dwarfs regime (a slope of -0.4 + 0.2 in logarithmic units). For 
both the upper Sco survey and for cr Ori, the mass function was 
determined over a large field of view, several square degrees and 
the associations are thus probed to much larger radii than e.g. 
previous studies of the ONC. However, although for some of the 
embedded clusters only the central regions have been observed 
(Mon R2, ONC, NGC 2024), others have much wider regions 
surveyed and the ratio of stars to brown dwarfs is still relatively 
high (e.g. Taurus, IC 348). 

The current finding of a possible change in the ratio of 
stars to brown dwarfs at larger radii suggests further studies of 
other nearby regions to larger radii than previously surveyed are 
needed. Also, obtaining kinematic data for the individual low- 
mass members to better characterise any mass segregation would 
be very valuable. Current near-infrared large field of view capa- 
bilities both from the ground and space are well suited for these 
studies. 

6. Conclusions 

We have presented deep NICMOS III F110W and F160W non- 
contiguous observations of the Orion Nebula Cluster region. The 
data reach 7=22 mag, and H- 21 mag, deeper than previous 
large scale high spatial resolution surveys. From the data we de- 
rive the ratio of low-mass stars (0.08-lM o ) to brown dwarfs as 
a function of radial distance from the centre. The ratio of stars to 
brown dwarfs is lower than in the centre, falling from a value of 
3.7 in the centre to 1.1+0.8 in the annulus 1.1-1.5 pc. This shows 
the low-mass content of the cluster is mass segregated to large 
radii. The age of the cluster is much shorter than the half mass re- 
laxation time suggesting the mass segregation for the low-mass 
population is primordial. The number of stars to brown dwarfs 
for the cluster as a whole is found to be ^ 03-008^) = 2.4 + 0.2, 
which is in better agreement with a flat or rising low-mass IMF 
than a Chabrier IMF. 
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Table 2. The ratio of stars to brown dwarfs in an extinction limited sample with Ay = Imag as a function of lower mass limit. Both 
the ratio corre c ted fo r fiel d stars and uncorrected is provided. For the inner 0.35 pc, the ratio of stars to brown dwarfs is taken from 
ISlesnick et alj d2004j) . The lSlesnick et al.1 ([2004) study only extends to 20 Mj up and the total ratio of stars to brown dwarfs down to 
10 Mj up is incomplete for the inner 0.35 pc. 



Corrected for field stars Uncorrected 

Ratio 0.35-0.7pc 0.7-1. lpc l.l-1.5pc Cluster 0.35-0.7pc 0.7-1. lpc l.l-1.5pc Field 

wff) ( 03 08 oo8) 7.2+5.6 1.6±0.8 1.1±0.8 2.4+0.2 f =4.9±1.8 § =1.6+0.5 g =1.4+0.4 f 

wffl02 8 oo8) 4.4±2.8 1.4±0.7 0.8±0.6 1.7+0.1 =2.9±0.9 § =1.2+0.3 § =0.9+0.2 ±§ 

^o (0 o| 8 o.o8) 4 - 1±3 - 3 i- 4 * - 8 1-2*1-2 1-6+0.1 ||=2.2+0.6 § =1.0+0.3 | =0.8+0.2 § 



